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ABSTRACT 

We use new large area far infrared maps ranging from 65-500 /im obtained with the AKARI and the 
Balloon-borne Large Aperture Submillimeter Telescope (BLAST) missions to characterize the dust 
emission toward the Cassiopeia A supernova remnant (SNR). Using the AKARI high resolution data 
wc find a new "tepirf" dust grain population at a temperature of ~ 35 K and with an estimated mass 
of 0.06 Mq. This component is confined to the central area of the SNR and may represent newly- 
formed dust in the unshocked supernova ejecta. While the mass of tepid dust that we measure is 
insufficient by itself to account for the dust observed at high redshift, it does constitute an additional 
dust population to contribute to those previously reported. We fit our maps at 65, 90, 140, 250, 
350, and 500 ^m to obtain maps of the column density and temperature of "cold" dust (near 16 K) 
distributed throughout the region. The large column density of cold dust associated with clouds seen 
in molecular emission extends continuously from the surrounding interstellar medium to project on 
the SNR, where the foreground component of the clouds is also detectable through optical. X-ray, and 
molecular extinction. At the resolution available here, there is no morphological signature to isolate 
any cold dust associated only with the SNR from this confusing interstellar emission. Our fit also 
recovers the previously detected "hot" dust in the remnant, with characteristic temperature 100 K. 
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1. INTRODUCTION 

Determining the origin of cosmic dust is fundamental 
to our understanding of many astronomical processes, 
including star formation and galaxy evolution. Galaxies 
and quasars at high redshift have b een found to contain 
large amounts of dus t (> 10^ M^ : IDunlop et al.l 119941: 
Archibald et~all 120011 : llsaak et al.l 120021: iPriddev et al.l 
20081) ■ at a time when the Universe was only about 
one tenth of its present age. The main source of 
dust injection within our Galaxy is thought to be 
the stellar winds of stars on the asymptotic giant 
branc h (AGB) of the Hertzsprun g-Russell (HR) dia- 
gram (jMorgan and Edmunds! [20031 ) . Stars at this early 
epoch would not have been able to reach the AGB 
phase in the available time, and therefore cannot be the 
source of the observed dust. Heavy elements are pro- 
duced in the explosions of supernovae (SNe) and, for 
several years, models have predicted that considerable 
amount s of fresh dust (0. 1 -1 Mf^) could also be pro- 
duced (Kozasa et a ll 119911: I Wooslev and Wc aved [l995l: 
iClavton et al. 199^ lTodinian d Fcrrara 2001). The life 
cycle of high mass stars (>8 M0), the progenitors of type 
II SNe, is sufficiently short for SNe to occur within the 
required timescales. As a result, SNe have been proposed 
as a possible solution for the origin of the dust seen at 
high- redshift. However, in order for SNe to generate suf- 
ficient dust mass to fill this ga p, each supernova would 
need to generate 0.4-1 M© dust pwek et al.ll2007[) . This 
quantification does not account for dust grain destruction 
within the supernova, thereby making it a lower bound. 

In seeking evidence regarding this hypothesis, the fo- 
cus has been on dust detectable in SNe and supernova 
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remnants (SNR), as reviewed briefly below. However, it 
seems less well apprecia ted that injection is only part of 
the story. iDraind (|2003[ ) reviews the often-ignored argu- 
ments that, at least in our Galaxy, the interstellar dust 
is continually processed on a timescale 3 x 10^ yr and 
most of its mass is formed in the interstellar medium. 
Nevertheless, injected dust is critical at least as seeds for 
further evolution of the dust population. 

Studies of SNe find only trace amounts of dust in 
the hot ejec t a, wit h typical ma sses of o rder 10~^M(n 



(iDwek et aU 119871: [Laeaee et al. 199& lArendt et all 
119991 : lErcolano et all 120071 : iMeikle eTaLl I2007D . Dust 



studies in SNR appear more promising, a prime ex- 
ample being Cassiopeia A (Cas A). Cas A is the rem- 
nant of a type li b supernova event which occurred 
around AD 1680 (iRavmondl 119841: iThorstensen et all 
[20Mt iFteSCTTeT ar '2006at iKrause et al.ll2008| ). The pro- 
genitor star is believed to have ha d a mass greater than 
20 Mq IE erez-Rendon et al] I2002D . and the remnant is 
at a distance D ~ 3.4 kpc (iReed et al.lll995D. Early ob- 
serva tions of Ca s A made with IRAS (jNeugebauer et all 
119841 ) and ISO (|Kessler et all 119961 ) did not extend to 
longer wavelengths, and therefore detected only the "hot" 
(^100 K) dust component, whose mass seemed insuf- 
ficient to provide the levels of dust seen in high red- 
shift galaxies. The low angular resolution also made 
the study of sub-structure in the remnant at interme- 
diate wavelen gths difficult. The Spitzer Spa ce Telescope 
(|Wcrnci''2006') has been us ed to study C as A (jHines et al.l 
[2004: .Krausc eTaP [200l lEnnis et al.l [20061 ). Most re- 
centlv lRho et al.l ( 2008() . exploiting the angular and spec- 
tral resolution achieved with the Spitzer infrared spec- 
trograph (IRS: iHouck et al.|[200l . find between 0.020 
- 0.054 M0 of hot dust. This is an order of magnitude 
greater than that previously measured and they conclude 
that, within modeling uncertainties for galaxy evolution, 
this could be sufficient to explain at least the lower limit 
to the dust levels in high-redshift galaxies presented by 
llsaak et all (|200l . 

Much more controversial is the question of a "cold" dust 
component in the SNR, because of the issue of contam- 
ination by linc-of-sight interst ellar emission from dust 
which is also cold (~16K; §0. iDunne et all (|2003[) . us- 
ing data from the Su bmillimetre Commo n User Bolome- 
ter Array (SCUBA: [Holland et al.l [19991) . find evidence 
for 2 — 4 Mq of dust in Cas A at a temperature of ~18 K, 
significantly more t han the mass of hot dust. Using the 
same mcthodologv. iMorgan et al.l (|2003[ ) find ~lMo of 
cold dust in Kepler's supernova remnant, a thousand 
times greater than previous measurements for this SNR. 
Both of these remnants are sufficiently young for the dust 
observed to be freshly formed in the remnant, rather than 
being material swept up from the I S M by the shock-wav e 
(jPickel et al.l U 
More recently. 



,9881: iHughes I [1991 IWright et all [19991 ). 
Dunne et al.l (|2009[ ) presented further ev- 



idence for cold dust in Cas A using SCUBA polarization 
data. These show dust emission polarized in an orienta- 
tion consistent with that of the magnetic field deduced 
from the radio synchrotron emission, suggesting the de- 
tected dust is in the SNR. They find a conservative lower 
limit for this dust mass of 1 Mq . They also attribute ap- 
parent depolarization at the brightest feature to dilution 
b y line o f sight interstellar emission. 
iDwekl ([2004i) has argued that the mass estimates in 



IDunne et al.l (|2003f ) exceed the total anticipated mass 
ejection of Cas A and suggest that if the subniillime- 
ter observations are valid they might instead imply the 
presence of a much smaller amount of dust which is a 
much more efficient millimeter wavelength radiator, such 
as iron needles. The alignment of such needles could af- 
fect the polarization. 

By correlating cold dust emission with molecular 
line absorption aga inst the SNR synchrotron emission, 
IKrause et al.l ()2004l ) argued that the dust is, in fact, as- 
sociated with a molecular cloud located along the line of 
sight to the SNR. Molecul ar emission has also bee n stud- 
ied in this direction (e.g., iLiszt and Lucaslll999D . show- 
ing that the c loud(s) extend well beyond the SNR itself. 
IKrause et al.l estimate the fresh dust yield within Cas A 
to be at least an order of m agnitude lower than that 
found bv IDunne et al.l (|2003f l. Even so, this would still 
provide the predominant dust mass in Cas A, bolstering 
the possibility of explaining the quantities of dust seen 
at high-redshift. 

The emission from cold dust peaks in the far-infrared 
and submillimeter and so neither SCUBA, on the long 
wavelength side, nor Spitzer, on the short side, is ideal 
for isolating a cold dust component. Near the ther- 
mal peak, the best large scale maps covering Cas A 
and its environs are from the Multiband Im aging Pho- 
tometer for Spitzer (MIPS: iRieke et al.ll2004l) at 160 /im 
(jKrause et al.ll2004[) and from the ISO Serendipity Sur- 
vey at 170/im (Sti ckel et al.ll2007f ). Neither map has full 
coverage, with MIPS missing data on small scales and 
ISOSS on larger. 

We observed Cas A us ing the Far-Infrared Surveyor 
(FIS: lKawada et al.|[2007[) instrument on-board AKARI. 
We obtained fully-sampled images of sub-arcminute res- 
olution covering a wide area surrounding Cas A in four 
photometric bands from 50 to 180 ^m. This is the wave- 
length range over which the emission from newly formed 
hot dust becomes faint and emission from cold dust be- 
gins to dominate. Therefore, these AKARI FIS images 
are very useful for investigating the presence of colder 
components of dust in the remnant. 

The Balloon-born e Large- Aperture Submillimeter 
TelescopeEB (BLAST; iPascale et all [20081 ) was also used 
to observe the Cas A region at 250, 350, and 500 ^m. 
These bands fill in the cold dust spectral energy distri- 
bution (SED) on the long wavelength side of the peak. 
The high mapping speed of BLAST means it was possible 
to cover a large area surrounding the SNR, giving data 
for both the SNR and the interstellar cloud structure in 
the surrounding region. 

By contrast, the prior SCUBA maps arc of a small area 
and involve deconvolution of a three-beam chopping pat- 
tern which reconstructs large scale power poorly. AKARI 
and BLAST have the advantage that, like Spitzer, they 
are not required to perform chopped observations. As 
a result, our new maps are sensitive to the large scale 
structure present in the Cas A field. The relatively large 
area of these maps, coupled with their sensitivity and 
wavelength coverage, make them ideal for investigation 
of cold dust emission from Cas A and the interstellar 
clouds. 

We describe the AKARI and BLAST observations in 

www. blastexperiment . info 
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§ m Several distinct sources of emission are distinguish- 
able in these maps. Taking advantage of the AKARI 
spatial resolution, in § |3] we identify a new morpholog- 
ically compact source of "tepid" (-^SSK) dust emission 
centered on the SNR whose spectrum is also distinct, 
peaking between the SED peak of the hot dust in the 
SNR and that of the cold dust. In §2] we perform aper- 
ture photometry on the maps in the region of the Gas A 
SNR. The resulting global SED further illustrates the dif- 
ferent spectral components and shows that without the 
additional morphological information it is not possible 
to unambiguously distinguish the tepid dust component. 
In § E] we fit the six-band AKARI-BLAST data with a 
simple spectral model to make column density and tem- 
perature maps for the cold dust. These clearly illustrate 
the confused nature of cold dust emission on the line of 
sight to the SNR. In i; 15. li the derived column density on 
the line of sight is compared to that obtained by other 
techniques. We present our conclusions in § [51 

2. DATA 

2.1. AKARI 

The AKARI FIS observations of Gas A were obtained 
on 2007 July 19 (Obs ID: 1402802). FIS has four pho- 
tometric bands at 65, 90, 140, and 160 /.im; two wide 
bands, 9 and 140 /im, and two narrow bands, 65 and 
160 /xm (jKawada et al.l [2007[) . FIS carried out a 'slow- 
scan' observation of Gas A with a scan speed of 15" s~^ 
and a reset interval of 1 s. The size of the detector pixels 
was 26'.'8 for the 65 and 90 fiui bands and 44'.'2 for the 
140 and 160 fim bands, approximately equal to the beam 
size. The area covered for Gas A is 40' x 12'. During the 
observation, the dark signal was measured five times, at 
the beginning, at the turning points of each scan leg, and 
at the end of the observation, by closing a cold shutter. 
The variation of responsivity was monitored by using an 
internal calibrator. 

The raw da ta were processed wi th the official FIS slow- 
scan toolkit (jVerdugo et al.l [20071 ). The removal of bad 
data as well as cosmic-ray hits and the correction of inte- 
gration ramps were made in the pipeline. During the re- 
sponsivity variation correction, dark current subtraction 
and flat fielding were performed. Since the dark level for 
160 /im was overestimated due to a large change of re- 
sponsivity, the 160 /im data were excluded in this paper. 
The sigma-clipping method was applied to calibrated sig- 
nals within the size of a detector pixel in order to remove 
small undetected glitches as well as other artifacts. The 
image was reconstructed with the assumption that a pixel 
value represents the uniform intensity over the pixel sur- 
face in the fine image grid. The average value was taken 
from the multiple values in the image grid. The abso- 
lute calibration was accomplished by comparison with 
data fror n the Diffuse Infrared Background Experiment 
(DIRBE; iSilverberg et al.lll993i ). For point-source cali- 
bration, the uncertainties of fiux calibration are no more 
than 20% for the two short wavelength b ands, and 30% 
for the 140 /im band (jKawada et al.ll2007f ). 

2.2. BLAST 

BLAST made a 'cap' observation (jPascale et al.|[2008l ) 
of Gas A during its first long duration balloon fiight 
from Kiruna, northern Sweden, in June 2005. Maps 



were made with a slew rate of 6' in azimuth, and 
a slow drift in elevation which produced a scan line spac- 
ing of 3f25 (about half an array width). Due to fiight 
time limitations, it was only possible to perform 1.25 
complete scan maps. The final images are ^^40' x 60' 
in size. The maps were made with the SANEPIG algo- 
rithm ([Patanchon et al.]|2008f) and were calibrated using 
the procedure discussed in iTruch "eFall ([2001 . 

Although the mirror diameter was 2 m, the BLAST05 
SANEPIG maps had only 3' full-width at half-power res- 
olution due to an anomalous telescope beam, corrupted 
by some uncharac terized combinati on of mirror distor- 
tion and de-focus (jTruch et al.ll2008l ) (this was corrected 
for the BLAST06 flight). Nevertheless, these maps have 
high signal-to-noise and are oversampled with 20" pix- 
els, so that Lucy-Richardson (L-R) deconvolution can be 
used to improve the resolution significantly (Roy et al. 
2009, in preparation). After 32 iterations, the L-R maps 
used here have point spread functions of lf3, 1.'6, and 
1'.9 in the 250, 350, and 500 /im bands, respectively, de- 
termined by measuring the size of a point source in the 
maps. While greatly improved, these are still not as small 
as the nominal BLAST values of 0.'5, 0.'66, and 1', or the 
even smaller values anticipated with the Herschel Space 
ObservatoriF^. 

The AKARI maps have a well-established zero point, 
but the BLAST SANEPIG processing fihers the lowest 
spatial frequencies and produces maps with zero mean 
value. In order to model the column densities and tem- 
peratures in § [5l we need to establish the true zero points 
of each BLAST map. As described below, we use DIRBE 
maps of the region to establish the appropriate mean 
value of the 250 /im BLAST map, and we use correlations 
between the BLAST maps to establish the zero points at 
the two longer BLAST wavelengths. In practice these 
two steps are done simultaneously. 

Although the DIRBE beams are 1.3 x lO^^'sr, much 
too large to resolve the structure around Gas A (there 
is of order one beam in the entire BLAST map), the ex- 
periment was explicitly designed to measure total flux 
and is therefore ideal for establishing the absolute inten- 
sity in each map. For the Gas A position the DIRBE 
intensities are 24.5, 86.8, 191.2 and 142.9 MJysr^^ at 
wavelengths of 60, 100, 140 and 240 /im. These values 
are plotted in Figure [T] 

The morphology of the structure in the BLAST maps, 
shown in Figure [21 are generally similar to each other, in- 
dicating that the SED of the interstellar dust is compara- 
tively uniform across the region. The Pearson correlation 
coefficients of the 350 and 500 /im maps to the 250 /im 
map are r = 0.98 and 0.97, respectively. The measured 
slopes. A, for the relation I\ = A/250 are 0.51 and 0.28, 
respectively, for A = 350 and 500 /tm, where BLAST map 
intensities, I\ have zero mean and have been measured 
on maps with the region of the Gas A SNR masked out. 
In the simple picture that the intensity variation within 
the maps is dominated by variations in column density 
from place to place, these linear relations should extrap- 
olate to the origin if the maps had correct zero points, 
and so knowledge of the absolute intensity in any one 

21 |http : // www . esa. int/SPECIALS/Herschel/index . html I 

http : //Icimbda . gsf c . nasa . gov/product/cobe/brows er . cf m| 
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Fig. 1. — Mean values of the maps within a DIRBE beam. The 
DIRBE photometry used to set the zero-points in the BLAST maps 
is shown as triangular data points, and the relative amplitudes of 
the BLAST mean values, as described in the text, are shown as 
diamonds. The solid curve is the sum of a cold dust curve with 
emissivity index /3 = 2 and adjustable temperature. Also shown is a 
curve for non-equilibrium dust emission from small grains with /3 = 
2, T = 40K and an amplitude 14% of that of the cold dust. The 
amplitude and temperature of the cold dust is chosen to minimize 
the deviation of this sum from the data points at the locations 
of the vertical error bars. The size of those error bars shows the 
assumed constant relative uncertainty corresponding to = 1 for 
the best fit. The pluses show mean values from the AKARI maps 
well away from the Cas A SNR, and the agreement is excellent. The 
inferred mean zero-point values appropriate to the three BLAST 
bands arc written next to each point. 

BLAST band would allow inference of the zero points of 
the other two maps. 

As described in § O however, there is a mild anti- 
correlation of temperature and column density: the dust 
in the dense arms of the molecular cloud is approximately 
2 K cooler than the dust in the less dense gas surround- 
ing them. This biases the measured slopes upward by 
7% and 11% at 350 and 500 /xm. 

The BLAST slopes, after correction for this small bias, 
are also plotted in Figure [1] relative to the 250 ^m point 
on the solid curve. This curve was obtained by fitting 
to the DIRBE fluxes at 100, 140, and 240 ^im and the 
relative BLAST amplitudes at 350 and 500 ^m, minimiz- 
ing while varying the amplitude of the curve and the 
temperature of the cold dust. We assume the same frac- 
tional uncertainty in each band. The solid error bars in 
the figure denote the bands used in the fit, and the size 
of the error bars has been adjusted post-fit to show the 
fractional uncertainty which would produce a reduced 
= 1. The cold dust temperature was found to be 
16.7 K. In order to accommodate the 60 fj,m data, the 
curve includes a contrib ution of non-equilibrium emis- 
sion ([Desert ct al. 199(|; ILi and D rainc 2001) from very 
small dust grains (VSGs) at an assumed temperature of 
40 K whose peak amplitude has been set to be 14% of 
that of the cold dust component. However, the 60 fim 
data are not used in the fitting procedure. Both the cold 
dust and the non-equilibrium components are assumed to 
have power law emissivities, e (x with /3 = 2. Chang- 
ing /3 to 1.5 does not improve the quality of fit, nor does it 
alter the BLAST zero points by more than a few percent. 

As a check, after the above fitting procedure was com- 
plete, the results were compared to mean values from the 
AKARI maps, which are shown as pluses in Figure [TJ 
The agreement is excellent. The BLAST maps presented 



in Figure [5] have had the correct zero points applied. 
2.3. Maps 

Figure [2] presents the AKARI 65, 90, and 140 /im and 
BLAST 250, 350, and 500 /im maps. All maps are cen- 
tered on the same location, arc smoothed to a common 
resolution of 1.'9, and are displayed with the identical 
range of brightness. 

The relative intensity of emission from the SNR and in- 
terstellar clouds varies significantly across the six bands. 
In the 65 /im image (Figure [2^) the SNR clearly domi- 
nates the map, with the cloud structure barely visible. 
Hot dust emission arises in the shell of the SNR where 
there is a reverse shock. This is evident in the shortest 
wavelength maps in the bright sources along a ring at 
the apparent SNR outer edge. To illustrate the spectral 
signature, the brightness at the location corresponding 
to the brightest peak in the 65 /im map minus the local 
background outside the SNR is plotted in Figure [31 A 
point at 24 /j,m is included to suggest the many short- 
wavelength data that are available (§111). When modeled 
as a simple single-temperature modified blackbody, the 
hot dust temperature Th is ~ 90 K (/? = 2). Note that 
this is a great simplification. Adopting a lower /? would 
broaden the spectrum, possibly accounting for an under- 
lying intrinsic temperature distribution. Nevertheless, 
Th is in the rang e found by the more detailed models by 
IRho et all (|200l . 

The SNR is also dominant in the AKARI 90 /.m image 
(Figurcl^b). However, the emission ceases to be confined 
to a simple ring. The brightness at the center of the 
SNR rises slightly going from 65 //m to 90 /im, which is 
not consistent spatially or spectrally with emission from 
either the hot or the cold dust component described be- 
low. In § [31 we present an analysis of the higher resolu- 
tion AKARI data to estimate the properties of this new 
source of intermediate temperature "tepid" dust emission 
(see resulting SED in Figure [31). 

The brightness of the interstellar emission increases 
at the longer wavelengths, demonstrating that the dust 
there is at a lower temperature than in the SNR. This 
"cold" dust emission peaks between the 140 and 250 /im 
data (Figs. [2t and[2l), and following a cold greybody 
SED recedes in the two longest wavelength maps (nev- 
ertheless, it is still detected with high signal to noise). 
The morphology of the cloud structure, both peaks and 
valleys, is clearest in the 140 - 350 /^m maps (Figs.[2fc - 
W^). Near the peak of the cold dust SED, the SNR does 
not stand out as a special feature in the map. There is 
a brighter cold cloud to the north-east in the extended 
BLAST maps. Closer to the SNR, there is a bright 
"ridge" of interstellar emission to the south-cast and an- 
other elliptical "cloud" beginning to the west outside the 
remnant and appearing to project across its face. These 
two features are labelled "R" and "C" in Figure [Hi re- 
spectively. It is clear that any emission associated with 
a possible component of fresh cold dust within the SNR 
will be highly confused by line of sight interstellar emis- 
sion; this particular interstellar cloud will be called the 
"line of sight cloud". The temperature of the cold dust 
throughout the map is ~ 16 K (§ [SJ. A representative 
SED is shown in Figure [31 for the molecular cloud ridge 
to the south-east of the SNR. 

The diffuse interstellar emission across the 65 - 90 fim. 
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Fig. 2.— AKARI 65, 90, and 140 Atm (a-c) and BLAST 250, 350, and 500 Atm (d-f) maps of the Cas A SNR and surrounding ISM, 
brought to a common resolution of lf9 and centered on the SNR. All maps are displayed with the same intensity scale. Contours are at 
the native resolution, in steps of lOMJysr"^. The letters C and R label the line of sight elliptical "cloud" and "ridge" features described 
in the text above. 



maps is generally too strong to be modeled along with the 
long wavelength emission as a single-temperature modi- 
fied blackbody. It is attributed to the above-mentioned 
non-equilibrium emission by VSGs 12. 2p . 

Synchrotron emission from relativistic electrons accel- 
erated in the SNR h as been mapped e xtensiv ely, includ- 
ing in the millimeter (jLiszt and Lucasl l999: Wright et al.l 
11999( 1. The synchrotron power-law SED rises smoothly 
toward longer wavelengths, and so eventually dominates 
over the falling dust spectrum. This is already becoming 
evident in our 500 ^m BLAST map, where the SNR has 
recovered in relative brightness compared to the inter- 



stellar emission. By 850 yum (SCUBA), the synchrotron 
emission is dominant. A representative SED is shown in 
Figure[3]for the synchrotron peak in the west of the SNR. 
This is drawn for a power-law s p ectral index a = —0.72, 
near that found bv iRho et al.l ()2003f ) . Spectral index 
variations across t he face of the SNR (IRho et alj 120031: 
Dunne et al.l l2009f l and with frequencv (jAtoyan et alJ 
2000( 1 are possible, but not critical for the analysis here. 



3. TEPID DUST INTERIOR TO THE SNR 

As discussed above, there is well-documented hot (~ 
100 K) dust in a broken ring closely associated with the 
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Fig. 3. — Illustrative SEDs of four distinguishable sources of emis- 
sion (not co-spatial). The "hot" dust spectrum, shown as triangles, 
is the intensity in each map taken at the location of the peak of the 
65 Mm map, o = 23''23'".4, <5 = 58°50'.l (J2000) minus the local 
background measured at 23''23'", -|-58°45'. The illustrative curve 
near these points is a /3 = 2 modified blackbody at T = 92 K. "Cold" 
dust (diamonds) is represented by the intensity of the "ridge" in the 
interstellar cloud structure to the south-east at 23''50'", -|-58°44', 
minus the respective map average. A T = 15K, /3 = 2 modified 
black body is drawn to guide the eye. The pluses show the in- 
tensity of a "tepid" dust Gaussian central contribution, determined 
in § [3] The two curves shown (/3 = 1, T = 41 K, dashed; and 
/3 = 2, T = 33 K, dot-dashed) fit the data equally well. The lower 
solid curve shows the synchrotron spectrum measured at 23''23™2, 
-|-58°48.' 8, the brightest peak in the west of the SNR in the 86 GHz 
map by ILiszt and Lucaa II1999I ) , and extrapolated using spectral 
index -0.72. 

rever se shock and the X-ray emitting plasma ()Rho et al.l 
|2008[ ). However, toward the core in the ejecta that has 
yet to be overrun by the reverse shock there could be 
cooler d ust that can be detected only at longer w ave- 
lengths pwek and Wernejll98lHMezger et allll986l ). Its 
minimum temperature would be set by the ambient ra- 
diation field and so it would be at least as warm as the 
surrounding interstellar dust. 

To track "tepid" or intermediate temperature dust 
which is neither the hot SNR ring component nor cold 
interstellar emission, the optimal spectral region is 90 - 
160 ^m. Indeed, using mul ti-wavelength (60 , 100, 170, 
200 ^m) ISOPHOT images. [Tilffs et al.l (|2005[) presented 
evidence for emission from such tepid dust in the un- 
shocked ejecta in two ways: (i) the morphology of the 
100 /im map shows more central emission than expected 
from the hot dust component which dominates at 60 /im, 
indicating a cooler centrally-peaked component, and (ii) 
the 170/200 brightness ratio is peaked on the SNR rather 
than extending with the band of interstellar emission well 
beyond the SNR to the west, indicating an emission com- 
ponent associated with the SNR that is warmer than the 
cold interstellar dust. We have confirmed these qualita- 
tive findings using the ISOPHOT pipeline-processed data 
from the ISO archive. 

The AKARI angular resolution and wavelength cov- 
erage are ideally suited for pursuing approach (i) to 
the next level, using a "spectrum-informed clean" tech- 
nique. Like the 60 /im ISO image, the 65 ^m AKARI 
image is dominated by the ring component. This can 
be isolated by subtracting the local background of about 
55MJysr~^. This "background-subtracted" map, a hot 
dust template, can then be scaled, following a hot dust 



SED, and subtracted from the background-subtracted 
90 /xm image. Choosing the SED so as to remove (i.e., 
"clean") the morphological signature of the hot dust ring, 
the striking result is an extended residual feature at 
90 ^m centered at 23'*23"29", -t-58°48'33" (J2000). To 
the accuracy available here, this direction is probably in- 
distinguishable from the direction to the center of the 
SNR. The shape of the feature is well approximated by 
a Gaussian of FWHM 2' which is substantially larger 
than the beam (37") though smaller than the SNR. The 
central surface brightness of this residual is about 20 
MJysr"^, which can be compared to the total of 133 
MJy sr~^ in this same direction. Following the same pro- 
cedure with the 60 and 100 /xm pair of ISO images pro- 
duces a residual of very similar size and position, with 
central surface brightness 27 MJy sr~^. This is an impor- 
tant check, because the images are not completely free of 
stripin g. This quantifies the emission component identi- 
fied bv lTuffs et al.l (|2005[ ) by its relative morphology. 

An interesting question is whether this extended Gaus- 
sian component can be seen in the short wavelength im- 
ages that have been used to generate the hot dust tem- 
plate. Depending on the three-dimensional geometry, 
there will be some hot dust emission projected on the 
central line of sight, though a peaked Gaussian arising 
in this way would seem unusual. Furthermore, at 60 and 
65 /im the surface brightness toward the center is not only 
brighter than the background, but compared to the hot 
dust ring is relatively brighter than is the case at shorter 
wavelengths such as in the archival MIPS 24 /im image 
(convolve d to the same reso lution) . This is also the con- 
clusion ofjRho^eta3j200l) based on the MIPS 70 /xm 
image (|Hines et al.ll2004[ )~~Therefore. to probe this more 
probable alternative, we assume that the central emis- 
sion in the background-subtracted 65 /im map can be 
attributed entirely to this Gaussian, which then would 
have a peak brightness of 33MJysr~^. 

Subtracting this Gaussian component from the 
background-subtracted 65 /im map produces a new hot 
dust template with a larger dip in the middle. There- 
fore, when this is subtracted from the 90 \m\ map, the 
Gaussian residual is still there, but with a larger ampli- 
tude, 42MJysr"^ (see Figure g]). Analysis of the 60 - 
100 /xm ISO pair yields a similar result. We consider this 
to be an upper hmit. The 60 - 65 to 90 - 100 /xm color 
temperature is about 33 K (/3 = 2) or 41 K (/3 = 1). 

Proceeding with this "spectrum-informed cleaning", 
now removing both the hot dust and tepid Gaussian com- 
ponents from the 90 or 100 /xm images, we find toward 
the SNR only a fainter residual in the west, (coinciden- 
tally) at the projected distance of the reverse shock, and 
extending to the southeast, a morphology crossing the 
remnant that is qualitatively consistent with being inter- 
stellar. The latter morphological component (called the 
"line of sight cloud" in § 12.31 and below) can be seen di- 
rectly in the 250 and 350 /im BLAST images (Figure d]), 
even without any component subtraction, and in the cold 
dust column density map (Figure [BJd below in § [5]) . As 
in the MIPS 160 /im image (pgrausc et al.ii2004l) and the 
long wavelength ISO maps, this extends to the west well 
outside the SNR. 

In the two intermediate-wavelength images (AKARI 
140 /xm in Figure[2}:; ISO 160 /im) the hot dust ring com- 
ponent is no longer obvious. But in both there are two 
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Right Ascension 
Fig. 4. — AKARI 90 background-subtracted map of the 
Cas A SNR from which the hot dust emission has been removed 
("cleaned"). Heavy contours from the AKARI 65 /im map show 
the position of the hot dust ring and the thin contours from the 
BLAST 250 /im map show the cold dust emission including the 
"line of sight cloud" and the "ridge" to the south-east. 

distinguishable peaks of emission, toward the center and 
the band from the western peak extending to the south- 
east. The Gaussian component could be as bright as 
35 MJysr"^ (with the second alternative hot dust tem- 
plate), fairly consistent with extrapolation of the tepid 
dust SED from shorter wavelengths. Gompared to the 
situation at 100 /im, the residual line of sight cloud pro- 
jected on the SNR and the adjacent "ridge" of interstellar 
emission to the south-east have become brighter, as ex- 
pected for cold dust. 

The peak amplitude of this central Gaussian feature, 
from the three maps where it can be distinguished, is 
plotted as pluses in Figure |31 Two modified blackbody 
curves plotted there (/? = 2, T = 33 K and (3 = 1, 
T = 41 K) describe the measured amplitudes equally 
well. Both indicate only a small flux from this com- 
ponent at longer wavelengths where the combination of 
coarser BLAST angular resolution and comparatively in- 
tense line of sight emission have not allowed direct mea- 
surement. 

As mentioned above, hot dust emission is correlated 
with X-ray line and continuum emission and infrared line 
emission from the shocked gas. It is therefore interesting 
to ask if there is any tracer of the plasma that is coin- 
cident with the central tepid dust emission. ISmith et al.l 
(|2009l) show maps of several lines from different ions trac- 
ing different physical conditions in the SNR. The [Si II] 
34.8 fim emission in their Figures 2 and 7 is from the 
unshocked ejecta in the interior of the SNR inside the 
reverse shock boundary (see their Figure 4). The peak 
position of this central emission appears to be consistent 
with the location of our Gaussian component. L. Rud- 
nick (private communication) finds from a map of the 
central low density areas convolved to the AKARI reso- 
lution that the FWHM is comparable too. 

Identification of this tepid dust Gaussian component 
clearly depends on clues from how the SNR morphology 
changes with wavelength. Using the higher-resolution 
multi- wavelength morphological information anticipated 
from Herschel, it will be very interesting to characterize 
this component in more detail for comparison with such 



(3-dimensional) tracers. 

3.1. Mass of Tepid Dust 

Dust emission maps record surface brightness lu which 
in turn depends on the dust mass column density Md 
through 

AU = I,/K,B,{Td), (1) 

where is the dust mass absorption coefficient, and 
Td is the temperature of the particular dust component. 
Since dust in this Gaussian component is assumed to be 
associated with the SNR, at distance D, integration over 
the face of the remnant results in 

Md = S,D^lK,B,{Td), (2) 

where is the measured fiux density and A^d is the 
derived dust mass. 

Adopting a central surface brightness 40 MJysr"^ at 
100 /^m, FWHM 2', and Td = 33 K, we get a dust mass 
0.055 X (30 cm^g~-^)/Kioo Mq. It would be a factor two 
lower for Td = 41 K. As is the case for the mass esti- 
mates for other dust components, the other large source 
of uncertainty comes from kioo- The normalizing value 
used here at 100 liva . is appropriate to interstellar dust 
()Li and Drainel[2b01[ ) but the actual value might reason- 
ably be different by a factor three. 

The derived mass for this new tepid component is com- 
parable to th at estimated pre viously for the hot dust 
component bv lRho et al.l (j2008D (see §[T|). Following their 
arguments and caveats, such dust yields could contribute 
to the dust masses seen in high redshift galaxies, but they 
are less than the required level 0.4-1 Mq estimated by 
iDwek et"al] ()2007f ) to account for the dust seen at high 
redshift. However, when combined with the previously 
reported hot and cold dust masses this makes supernova 
remnants a more plausible source of this dust. 

4. GLOBAL SPECTRAL ENERGY DISTRIBUTION 

We measure the spatially-integrated ("global") flux 
density Si, in the direction of Gas A using a simple 
aperture photometry method. A three-arcminute radius 
aperture is used to measure the SNR flux, and an annu- 
lus with inner and outer radii of 4' and 5' respectively 
is used to measure the background. Only information 
from the western side of the annulus is used to estimate 
the background, as the eastern side contains elements 
of large scale cloud structure, which would result in an 
overestimate of the background flux. 

In addition to measuring the fiux density of Gas A in 
the six bands 65, 90, 140, 250, 350, and 500 ^im from 
AKARI and BLAS T, we remeasured t he Infrared Space 
Observatory (ISO: IKessler et al.l[l996h 170 /im photom- 
etry using data from the archive. We obtained a value 
of 101±2 0Jy, significantly higher than the 35±10 Jy ob- 
tained bv I Tuffs et al.l (|1999f ) which we believe was a lower 
hmit. The AKARI, BLAST, and ISO values are given in 
Table 1, and plotted in Fig ure [5] along with p ublished 
MSX, IRAS and MIPS data (jHines et al.ll2004l ). 

In order to fit these data we developed a physically 
motivated model with four dust components as discussed 
above, including hot dust (subscript H), the tepid Gaus- 
sian component (G), cold dust (C), and VSGs (V), to- 
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Flux Density (Jy) 


65 


71 ± 20 


90 


105 ± 21 


140 


92 ± 18 


170 


101 ± 20 


250 


76 ± 16 


350 


49 ± 10 


500 


42 ± 8 



TABLE 1 

Table of Cas A flux density measurements for the AKARI 

and blast bands, and the disputed iso 170 /im band. 
These were made using the aperture photometry described 
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Fig. 5.— Global SED of the Cas A SNR. A fit was made to 
publi shed MSX, IRAS, and Spitzer/MIPS photometry l IHines et aTl 
120041 1 together with the new AKARI and BLAST and remeasured 
ISO 1 70 fim (open circle) data, but neither the original lTuffs et al.l 
(119991 ) value (filled circle) nor the SCUBA data (triangles) were 
used. The model used for this fit is given in equation |3] and con- 
tains physically motivated components attributable to line of sight 
interstellar dust, and dust and synchrotron emission from the SNR. 
The two broad peaks in panel (a) are from hot SNR dust (dotted 
curve) and cold dust (upper dashed curve) with derived temper- 
atures 99±3K and 16.6±1.3K, respectively. More minor compo- 
nents are the new tepid dust Gaussian component in the SNR (§|3] 
dot-dash curve) and VSGs (lower dashed curve). These contribute 
in the spectrally-confused 100 /^m region and cannot be separated 
using this global SED alone. For panel (b), the flux of the syn- 
chrotron and hot and tepid dust components have been subtracted 
to highlight the contribution of cold dust and, less importantly, 
VSGs. 



gethcr with synchrotron emission, and described by: 

S.^Yl ^^^^'MTO + So{v/voT, (3) 



where the scale factors Q, are related to the spatially- 
integrated column density or mass. All fom" dust com- 
ponents have been approximated as single-temperature 
modified blackbodies. We obtain t he synchrotron fluxes 
by extrapolation from the map of iWright et al.l ()1999l ) 
and subtract them from Si, and so these data are not 
part of the fit to find the dust parameters. 

For the cold dust ( largely interstellar), we take f3c — 2 
(|Li and Drainel[2i301[ ) and find the other two parameters, 

and T, from the fit. 

The value of /3h is much less certain. Furthermore, 
this is a highly simplistic model, as the hot modified- 
blackbody underestimates the emission from polycyclic 
aromatic hydrocarbo ns (PAH) and short-wavelength 
SNR features seen bv iRho et all (|2008[) . /3h = 1 might 
be a better approximation to broaden the model SED to 
accommodate the range of hot dust temperatures, but we 
take it to be 2 to better describe the longer wavelength 
emission near 100 /im (Figure [3|). The hot dust emission 
from the SNR dominates the total flux emitted in the 
short wavelength region of the SED and so its character- 
istic temperature Th is well determined, but discovering 
exactly how the SED extrapolates to longer wavelengths 
requires morphological information, and is not the focus 
here. 

The single-temperature approximation is least appro- 
priate for the non-equilibrium emission by VSGs which 
would have a range of temperatures, but this is not 
critical here since the VSG emission is relatively weak. 
To help constrain the model, we take a representative 
Tv = 40 K, /3v = 2, and Hy/f^c = 1.8 x IQ-^ determined 
from the large scale interstellar regions of the maps (§[5]). 
Consistent with the integrated spectrum shown in Fig- 
ure [U this results in the peak of the VSG emission being 
about 14% of that of the cold dust. 

We take the parameters of the Gaussian component 
from our study in § [3l where the flux at 100 /im, near 
the peak, is found to be 16 Jy and the SED is as shown 
in Figure [31 

Fitting the remaining parameters to the global spec- 
trum, the total model spectrum and its components are 
plotted in Figure [5l The original ISO 170 /.jm point and 
the SCUBA data were not included in our flts. The ISO 
flux we measure is consistent with the new BLAST and 
AKARI data, as well as with previously published data 
at other wavelengths. This fit does not appear consistent 
with the SCUBA data at either 450 or 850 ^m, with the 
BLAST 500 Mm point being -40% (1.8cr) lower than the 
SCUBA 450 /im value. Wc cannot find any explanation 
for this disparity. 

We find Th = 99±3K and Tc = 16.6±1.3K. For the 
same k as adopted for the tepid dust (for lack of bet- 
ter constraints) , the independent masses of hot and cold 
dust (if the latter were all at 3.4 kpc) are 0.003 and 
10 Mq, respectively. The latter is the potential amount 
of cold dust associated with the SNR, significant back- 
ground emission having been subtracted when carrying 
out the photometry. This limit, close to that estimated 
bv iKrause ct al., (,2004i ) . illustrates just how difficult it 
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is to separate out eold dust in the SNR, given a strong 
interstellar signal. This too needs to be addressed using 
morphologieal information. It is noting that when the 
newly measured ISO point is exclulded from the fit there 
is only a negligible effect, changing the derived param- 
eters by less than l-tr as compared to a fit performed 
without this datum. 

5. COMPONENT SEPARATION 

The photometry in Figure[5]benefits from greater spec- 
tral coverage to provide an improved measurement of the 
spectral energy distribution toward Gas A, but does not 
explicitly address whether that emission, particularly the 
cold dust, is line of sight interstellar emission or is in- 
trinsic to the SNR. Here we exploit the six large area 
AKARI and BLAST maps over the region observed in 
common. The maps were convolved to a resolution of 
1.'9 and re-gridded to the same pixelization to produce a 
coarse spectral "cube". Using a multi-component model 
akin to equation [3] in §|11 an independent fit to the SED 
at each pixel in this cube is carried out and maps of each 
free parameter obtained. 

We first subtracted the SNR synchrotron component 
from the data and so, as in § 01 this model component 
is not a part of the fi t. We convolved the 83 GHz radio 
map of IWright et al.l (|19990 to produce a synchrotron 
template, and subtracted the appropriate intensity by 
extrapolating the spectrum assuming a power-law spec- 
tral index of a = —0.72 (see ij 12. 3p . This left us with 
maps containing only thermal dust emission, but still 
four components on the line of sight to the SNR, which 
introduces too many free parameters. Therefore, bene- 
fiting from the higher resolution analysis in § [3l we also 
subtract from the maps the tepid dust Gaussian compo- 
nent found there. 

Because the analysis here uses only the AKARI and 
BLAST maps, it does not contain information from wave- 
lengths shorter than 65 /im. Gonsequently it is not neces- 
sary or possible to have a sophisticated modeling of short 
wavelength emission sources (e.g., PAHs and shorter 
wavelength VSG emission). We simply used Ty = 40 K 
and /3v = 2 as above. Additionally, toward the SNR 
the intensity of the VSG component used in this simple 
model is degenerate with the hot dust ring component 
and tepid dust Gaussian component because of the over- 
lapping SEDs, meaning that our analysis is not sensitive 
to whether the ratio of VSG to cold dust emission for pu- 
tative SNR dust might be different than the equivalent 
ratio for the interstellar line of sight cloud. Away from 
the SNR we allowed this ratio to vary in order to ex- 
amine its statistical properties, and then we adopted the 
average ratio for the confused line of sight to the SNR. 
The properties of the cold dust (column density and tem- 
perature) outside the SNR are quite similar whether this 
average ratio is used in the model for the whole map, or 
is allowed to be fit. 

The hot dust component should only be included in the 
model for pixels where the SNR contributes to the inten- 
sity. In order to assess in an objective automated way 
whether or not the data to which the model is fit con- 
tain emission from the SNR, the AKARI 65 to 140 /im 
ratio was examined. A maximum ratio of 0.22 was found 
for all regions beyond a 4' radius from the SNR. There- 
fore, only when the ratio exceeded this value was this hot 



component included in the model fit. The temperature 
of the hot component is held fixed at 100 K, as derived 
in §[4] with /3h = 2, and its amplitude (column density) 
allowed to vary. As an alternative, we experimented with 
removing this emission component from the map using a 
scaled hot dust template. This final cleaning step was not 
entirely satisfactory, probably because of changes in the 
hot dust SED (temperature) across the face of the SNR, 
but again the derived properties of the cold dust were 
not very sensitive to this choice. As another sensitivity 
test, we fit using a hot dust model with Th = 120 K and 
/3h = 1 and again the cold dust parameters were robust, 
because the overlap of the cold and hot SEDs is not great 
(Figure p. 

Our six-wavelength data are fit using a model which 
contains either two or three free parameters, depending 
on whether the data being fit contains a contribution 
from the SNR. These parameters are the hot dust am- 
plitude, the cold dust amplitude, and the cold dust tem- 
perature, Tq. Implicit in this is the assumption that the 
cold dust model component is sufficient to describe the 
emission from both the SNR and line of sight interstellar 
clouds. A reduced measurement was calculated for 
each pixel fit; the median value over the map was 0.89. 

The output maps are presented in Figur e El with con- 
tours from the radio synchrotron emission (jWright et al.l 
Il999f ) overlaid to show the location and approximate size 
of the SNR. 

We find interstellar cloud temperatures ranging be- 
tween 15 and 18 K, with a median of 16.5 K (see Fig- 
ure [S^), in keepin g with the temperature found by 
iWilson et al.l (I1993D an d the cold dust temperature found 
bv lDunne et al.l (|2003D . The fine structure in this tem- 
perature map reflects some interplay between the various 
model components contributing to the emission model, 
but the main features there are robust against exactly 
how these spectral components are treated. 

The Gas A SNR does not stand out particularly. There 
is a slightly hotter spot to the north of the SNR whose 
amplitude changes with the choice of alternative hot dust 
models and so probably reflects a deflciency there where 
the hot dust component is very bright. There is no tell- 
tale change in temperature at the location of Gas A, to 
indicate the presence of a significant secondary cold dust 
grain population at a different temperature to the inter- 
stellar cloud. 

The cold component column density map (Figure ^jp) 
shows distinctly the cold dust ridge (R) to the south-east 
of the SNR, and the cold elliptical region (G) to the west 
of and overlapping the SNR (the "line of sight cloud") . 

Note that feature G is not as strong as R, is not obvi- 
ously correlated with any morphological features of the 
SNR. and extends to the west beyond the remnant. The 
cold column density map resembles what is seen in in- 
tegrated moleoila^Jin^emis^^ (Figure [71 see also Fig- 
ure 4 of iLiszt and Lu caslll9 99|. Figure 1 of iKrause et al.l 
[200l Figure 1 of [Wi lson a nd Batrlall2005l and Figure 8 
oflo unne et al.ll2009[ ). Therefore we conclude that there 
is, by chance, a substantial amount of interstellar ma- 
terial projected on the Gas A line of sight. This makes 
it difficult to make a statistically significant detection of 
cold dust distinctly associated with the SNR with these 
data. It is likely that this conf u sion a ffected the measure- 
ments made by [Dunne et al.l ([20031) based on chopped 
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Fig. 6. — Output maps for the three free parameters of the spectral model: temperat ure of the cold dust and column density of both cold 
and hot dust. In the upper panels a single white contour from the 83 GHz radio map llWright et al.in"999l ) shows the approximate position 
of the Cas A SNR forward shock IjHelder and ViiiHl2008l ) . The letters C and R label the elliptical line of sight "cloud" and "ridge" featured 
high-lighted in[2jl respectively. The contour lines vary in steps of 0.5K, 0.5 X lO^^cm"^, and lxlO~*gcm~^ for maps a-c, respectively. 
Note that the outer edge regions are corrupted by convolution artifacts. The lower panels show magnified versions of the output within 
the white dashed square centered on 2.3''23'"24^ -f-58°48'54", overlaid with several white radio contours. 

to the north west which intersects with the remnant. 
Comparison with the cohimn density map shows this to 
be part of a general anti-correlation between temperature 
and column density across the map. This could reflect a 
deficiency in the modeling, since the two are certainly in- 
versely coupled, non-linear ly, through equation [TJ But it 
seems plausible that this has a physical origin. The high 
column density regions are molecular and, being more 
shielded from the interstellar radiation field which heats 
the dust, would be cooler. The low column density re- 
gions also tend to have a higher ratio of VSG emission, 
when this is allowed to vary in the fit, though this might 
be indicative of another effect as well, evolution of the 
grain size distribution. 

The third output map. Figure [St, shows the column 
density from the hot SNR dust component. The auto- 
mated fitting has invoked this component only at the 
location of the SNR, the interstellar dust model being 
sufficient at locations away from the SNR. As expected, 
the structure of the emission is similar to that in the 
65 fim band when convolved to the same resolution as 
the BLAST data. 

Quantitatively, the hot and cold component ampli- 
tudes were converted to dust mass column density us- 
ing equation [1] In principle different values of k^, would 
be used for the different dust components in differ- 
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Fig. 7.— Integrated ^^C O (1-0) Une emiss ion from the FCRAO 
spectral cube in the CGPS l ITavlor et al.l2003 ). Thin contours from 
the BLAST 250 fira map show the cold dust emission including 
the "line of sight cloud" and the "ridge" to the south-east. Thick 
contour indicates the position of the radio SNR as in Figure (6] 

SCUBA observation s. We can not ident ify the polarized 
dust source found bv lDunne et al.l ()2009D : however, based 
on our analysis of the data, this is to be expected as the 
flux density level predicted would be too faint to detect 
in our observations. 
There is a large warmer-than-average region of '^17 K 
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ent environments but, in the absence of persuasive ev- 
idence, we proceeded as in §§[3]|4]and adopted those of 
iLi and Drai^ (l200l (kiqo = SOcm^g-i; /3 = 2) which 
are suitable for the diffuse ISM. Since much of the ISM 
in this region is in molecular form, these might not be 
entirely appropriate even for that component. The value 
of Ki, represents the greatest uncertainty in this analy- 
sis, with published values spanning two orders of mag- 
nitude, depending on the dust type and environment 
(|Dunne et al.l[2003[ ). 

For the cold dust component, which is dominated by 
the line of sight interstellar dust, we further convert 
the dust mass colu mn density into -/Vh by dividing by 
1.9 X 10^26 gjj-i ()Li and Draind I2001D . iVn, used in 
Figure |6}d, is the metric quantifying X-ray photoelectric 
absorption, discussed in § 15.11 The column densities are 
large throughout the map, of order 10^^ cm^^. Simi- 
larly, this can in turn be converted to opti cal extinction 
Ay b y dividing by 1.9 x 10^^ cm~^ mag~^ (jBohlin et al.l 
|1978() . equivalent to multiplying the mass column den- 
sity directly by 2.8 x 10^ cm^ g""'^ mag. The relationship 
between iVu and Ay is empirically cali brated only below 
about Nu = 0.4 x 10^^ cm^^ or Ay = 2 (|Kim and Martini 
|1996[) but might still be a reasonable approximation for 
the interstellar material spread out along this long line 
of sight. 

Note that both A^h and Ay, being scaled from es- 
timates based on optically-thin submillimeter emission, 
measure the total column extending through the Galaxy. 
The column in the foreground of Gas A should be compa- 
rable to this, simply given the large distance to the SNR 
and Galactic latitude, —2.1°, and longitude, 111.7°. In 
fact, the velocity range of molecular lines seen in ab- 
sorption against Gas A is the same as that seen in emis- 
sion (§ 15. ip suggesting that Gas A is beyond most of the 
Perseus arm gas. Judging from the H I 21-cm emission- 
line spectru m along adjacent l ines of sight in the CGPS 
data cube (jTavlor et al.ll2003| ). up to 0.1 x lO^^ cm^^ 
might be in the background, or less than 10%. 

If the cold dust projected on the remnant were all 
at that distance, the mass column density integrated 
over its face in a 5' diameter aperture would amount 
to ^ 40 Mq. This is a few times larger than found 
in § H] because here there has been no subtraction of a 
"background"; Figure [6b shows that the column density 
on surrounding inte rstellar lines of sight is substantial. 
iDunne et al.l (I2009D sugges t detection of 1 Af© in the 
SNR whereas iKTause et al.l (2004) find an upper limit of 
0.2 Mq. One has to beware of the different k and Tc that 
have been adopted in different derivations (the latter was 
not well constrained previously). In any case, it is again 
clear that the line of sight interstellar emission is an over- 
whelming source of contamination, with the searchcd-for 
SNR signal at the few percent level. Identifying cold SNR 
dust will be challenging even with the improved resolu- 
tion of Herschel working close to the peak of the cold 
dust emission. 

5.1. Other Measures of Interstellar Column Density 

A supplementary approach would be to look at other 
measures of the ISM column density as a surrogate for 
the morphology and brightness level expected for the con- 
taminating interstellar dust emission. 

One potential surrogate for column density is molec- 



ular line emission. Molecular emission is affected vari- 
ously by optical depth, abundance, and excitation con- 
ditions, and so a precise one-to-one correspondence with 
dust emission would not be expected. Nevertheless, even 
Figure |7l the integrated line emission from the FGRAO 
^fCO (1 - 0) data cube in the GGPS (jTavlor et al.l 
|2003() . looks remarkably similar, though not identical, 
to the cold dust emission over this extended region (Fig- 
ures [2]:-f ) and to the column density map (Figures [6j) 
an d e). Near the SNR, th e integrated GO (2 — 1) map 
by iLiszt and Lucad (|1999( ) also shows features recogniz- 
able in the dust emission, including the western peak 
and th e separated south-east ridge. iWilson and Batrlal 
()2005D emphasize how ^^CO (1 — 0) also projects on the 
SNR. These tracers suggest, as does the larger GGPS 
map, that the dust emission from the interstellar cloud 
G should extend beyond the SNR to the west. Note also 
the molecular emission toward the center of the SNR in 
these maps. The correlation of dust emission and molec- 
ular emission is not and cannot be expected to be perfect, 
but it is clear that interstellar line of sight dust emission 
is a major contaminant. 

If the interstellar material seen along the line of 
sight is substantially in the foreground of the Gas A 
SNR emission, as appears to be the case, then mea- 
sures of absorption can be used to trace column 
density. Ther e are three approaches: (i) molecu- 
lar lines (OH: [Bieeing and Crutcherl 119861 : methanol: 
Revnoso an d Goss' '2002*) and the H I 21-c m line (e.g., 
Keohane et aLniggftl . (ii) optical re ddening (iFesen et al.l 
2006bD. and (iii) X-ray absorption (iKeohane et al.lll998l: 
Vink et all [l99a IHuehes et all 120001: IWillingale et al.l 
20021) . These give qualitatively the same results: (i) the 
column density is highest toward the west of the SNR, 
where we see the residual at 100 and 140 /im after sub- 
tracting the tepid Gaussian component (§[3|), (ii) there 
are patchy enhancements along a band extending from 
the west to the south-east, thus covering the south-west 
and south portion of the remnant, like feature G in our 
cold column density map and the molecular emission, 

(iii) the absorption is lowest across the north-east, but 

(iv) there is a substantial column everywhere, including 
toward the center (and the neutron star), and (v) the 
surrogate column density is about the same as derived 
in the above component separation, which indicates that 
our choice of kioo was reasonable. 

Quantitatively, column densities from the extinction 
surrogates are of order 10^^ cm~^, in these units about 
2.3 toward the western maximum and 1.3 toward the 
center {Ay ^ 7)E3. Not only are these values and the 
low contrast in column density comparable to what we 
found independently above, but there is also morpho- 
logical similari t y. Go mparc, for exam ple. Figure Ic in 
IKeohane et al.l (|1996( ) and Figure 4 in IWillingale et al.l 
()2002|) to our cold dust column density map in FigurelHt. 
This leaves little room for an additional contribution at 
long wavelengths from the Gaussian "tepid" SNR dust 
componen t or a colde r one. 

[Krause et al.l ()2004[ ) show that the patchy OH absorp- 

Note that these column densities arc large enough to have 
an extinction effect on the infrared spectra of the SNR, although 
no absorp t ion (e .g., from silicates) was detected or modeled by 
IRho et al.l pOOSh . 
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tion correlates well with the residual 850 /im emission 
seen with SCUBA (after subtracting the synchrotron 
emission), and on subtracting this foreground emission 
(scaled using Td — 14 K) from their 160 /im map they 
find only a small residual and hence the above-mentioned 
upper limit of 0.2 M© of cold dust associated with the 
SNR. 

Accounting for and removal of the (foreground) emis- 
sion by correlations with extinction or molecular emission 
does seem a promising technique. It will be interesting 
to see if it proves useful pixel by pixel with forthcoming 
high signal to noise data or only statistically. Clearly the 
higher resolution imaging anticipated with Herschel will 
be essential for making progress along these lines. 

6. CONCLUSIONS 

1. We presented far-infrared/submillimeter data at 65 
- 500 /xm for the Cas A supernova remnant and the 
surrounding region. We used these maps to charac- 
terize the interstellar dust emission using data from 
cloud regions well beyond the SNR. 

2. We used high resolution ARAKI data to probe the 
spectral region between the hot dust emission from 
the SNR shock-front and cold dust emission. Us- 
ing a spectrum-informed clean technique we identi- 
fied a new tepid dust population at temperature of 
^35 K. The mass of this individual dust population 
was estimated to be 0.06 M0, but with considerable 
uncertainty because of its dependence on the choice 
of K. 



tentially significant source of dust production in the 
high-redshift universe. 

4. We developed a simple physically-motivated model 
of the SED of the SNR and interstellar emission 
and fit this to six-wavelength bands at each pixel. 
From this we obtained temperature and dust mass 
column density maps. The interstellar dust was 
found to be at a temperature of 16.5 K, in keep- 
ing with previous measurements, but now better 
constrained due to the improved wavelength cover- 
age. 

5. We show that the high level of confusion arising 
from the interstellar cloud structure projected on 
the SNR precludes a significant detection of cold 
dust directly associated with Cas A. The same 
source of confusion will have affected previous es- 
timates of cold dust in Cas A, increasing the un- 
certainty of those estimates. This analysis was not 
sufficiently sens itive to identi f y the lower limiting 
mass fou nd bv iDunne et al.l (j2009() or the lower 
value by iKrause et al.l ( 2004) and therefore does 
not preclude the possibility of a significant popu- 
lation of cold SNR dust grains with temperature 
close to that of the interstellar dust. The higher 
angular resolution data anticipated with Herschel 
working close to the peak of the cold dust emis- 
sion, together with correlations with surrogates of 
the interstellar column density, could result in a 
more sensitive probe. 



3. The dust yield for this new and independent tepid 
component is comparable to that esti mated pre- 
viousl y for the hot dust component by iRho et all 
(20081). While such yields could contribute to the 
dust masses seen in high redshift galaxies, they are 
stil l less than the requ ired level 0.4-1 MQCstimated 
by iDwek et al.l ([2007h . While the mass we mea- 
sure is insufficient to account for the dust observed 
at high redshift, when taken in combination with 
the hot and c o ld du st ma sses previousl y repo rted 
by IRho et al.l (|200l and IDunne et al.1 (pOOl ). it 
strengthens the argument for supcrnovae as a po- 
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